Pinned modes in lossy lattices with local gain and nonlinearity 
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We introduce a discrete linear lossy system with an embedded "hot spot" (HS), i.e., a site carrying 
linear gain and complex cubic nonlinearity. The system can be used to model an array of optical 
or plasmonic waveguides, where selective excitation of particular cores is possible. Localized modes 
pinned to the HS are constructed in an implicit analytical form, and their stability is investigated 
numerically. Stability regions for the modes are obtained in the parameter space of the linear gain 
and cubic gain/loss. An essential result is that the interaction of the unsaturated cubic gain and 
self-defocusing nonlinearity can produce stable modes, although they may be destabilized by finite- 
amplitude perturbations. On the other hand, the interplay of the cubic loss and self-defocusing 
gives rise to a bistability. 



I. INTRODUCTION 

Dissipative spatial solitons, which originate from the 
interaction of diffraction, self-focusing nonlinearity, and 
dissipation (gain and loss), have drawn great interest 
in nonlinear optics pQ and, more recently, in plasmon- 
ics [2j [3]. A necessary condition required for the for- 
mation of stable solitons is the stability of the zero 
background. The simplest single-component complex 
Ginzburg-Landau (CGL) equation with the uniform lin- 
ear gain is not an appropriate candidate for modeling dis- 
sipative solitons, as it violates this condition. Dissipative 
solitons may be stabilized by a system of linearly coupled 
CGL equations [4] that models dual-core waveguides with 
the linear gain and loss acting in different cores [3] HH7] . 
Stable solitons can also be generated by the single CGL 
equations that incorporate the linear loss, cubic gain and 
quint ic loss, which accounts for the nonlinear saturable 
absorption [HJ[9]. In these models, the quintic loss satu- 
rates the growth induced by the cubic gain and therefore 
stabilizes the solitons. 

Another method for generating stable localized modes, 
which has recently drawn considerable attention, relies on 
the action of linear gain at a "hot spot" (HS, i.e., a local- 
ized region in a lossy waveguide [TQUT2] or in a dissipative 
Bragg grating [13]). Models with several hot spots [14]- 
fT6] . as well as with similar extended structures [17 , have 
also been studied. HSs can be created by implanting an 
appropriate distribution of gain-producing dopants into 
the waveguide [18] , or by illuminating a uniformly-doped 
waveguide with external pump beam(s) focused at the 
designated spot(s). Dissipative solitons pinned to HSs 
can be stabilized via the balance between the local am- 
plification and uniform loss in the waveguide. In partic- 
ular, solutions for dissipative solitons pinned to narrow 
HSs approximated by delta-functions have been found 



analytically [TUl [14] [16]. Other relevant modes, both 
one- and two-dimensional, including stable vortices sup- 
ported by the gain applied to a confined area [19], have 
been generated in the numerical form [TTJ [12] [15] . Stable 
dissipative solitons have also been predicted in a system 
that combines the uniformly-distributed linear gain and 
nonlinear loss growing towards the periphery faster than 
r D , where r and D represent the distance from the center 
and spatial dimension, respectively [20] . 

An interesting ramification of the configurations men- 
tioned above is the possibility to generate stable solitons 
supported by localized cubic gain, in the absence of the 
quintic gain saturation. While dissipative solitons cannot 
be stable without higher-order nonlinear losses in uni- 
form media [2TJ [22], it was recently demonstrated [12] 
that stable dissipative solitons may be pinned to an HS 
carrying the unsaturated cubic gain and embedded into 
a uniform linear-loss background. This finding suggests 
ways to design clean nonlinear soliton amplification that 
avoids concomitant generation of noise, which is also rel- 
evant for plasmonics [23] . 

The present work explores the generation of stable 
solitons in discrete waveguiding arrays (lattices) with a 
localized unsaturated nonlinear gain. In particular, we 
demonstrate that this is possible in a linear lattice where 
the nonlinearity, represented by the self-phase modula- 
tion and cubic gain, is applied to a single waveguide (the 
HS). The lattice CGL system is introduced here as a dis- 
crete counterpart of the continuous HS models [10] [12] , 
and can be used to investigate various effects in pho- 
tonics, cf. Refs. [24] [27] [28]. In particular, it may be 
used for selective excitation of particular core(s) in an 
arrayed waveguiding system, if it is uniformly doped, but 
only the selected core is pumped by an external coherent 
source of light. In addition to allowing the straightfor- 
ward experimental implementation [24], we demonstrate 
that the present discrete system supports analytical solu- 
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tions for discrete solitons (similar to those in the discrete 
linear Schrodinger equation with embedded nonlinear el- 
ements [29]), thus providing a natural platform for ex- 
ploring the soliton dynamics. 

The paper is organized as follows. The discrete CGL 
equation with HS, and its implicit analytical solutions for 
pinned modes, are introduced in Sec. [TTJ The linear sta- 
bility analysis of the solitons against small perturbations 
is presented in Sec. |III[ and results of numerical computa- 
tions of the corresponding eigenvalue spectra are reported 



in Sec. IV The predictions of the linear-stability analysis 
are compared to direct simulations of perturbed evolution 
of the discrete solitons. In particular, stable solitons are 
found under the unsaturated nonlinear local gain, pro- 
vided that the localized nonlinearity is self-defocusing, 
although finite-amplitude perturbations may destabilize 
these modes. On the other hand, the interplay of the 
self-defocusing nonlinearity and cubic loss gives rise to 
bistability of the pinned modes, which is a sufficiently 
interesting finding too. The paper is concluded by Sec. 

El 



II. THE MODEL 

The present work is motivated by the complex 
Ginzburg-Landau (CGL) equation that models the prop- 
agation of an electromagnetic field of amplitude u(x, z) 
in a lossy waveguide with an embedded HS: 



du 

dz 



i d 2 u 
2£h? 



■ + [(r x + iY 2 ) + (iB - E) \u\ 2 } S(x)u . 

(1) 

Here x and z are the transverse and longitudinal coor- 
dinates, and 7 > is the background linear loss. Delta- 
function 5(x) approximates the concentration of the lin- 
ear gain (Ti > 0), linear potential (T 2 > corresponds 
to the local attraction) , cubic dissipation E (positive and 
negative for the loss and gain, respectively), and Kerr 
nonlinearity B (positive and negative for the self-focusing 
and defocusing) at the HS. 

As the model of an array of guiding cores, replacing a 
single continual waveguide, Eq. ([!]) is substituted by its 
discrete version, 
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where m = 0, ±1, ±2, ... is the discrete coordinate, 5 m ^ 
is the Kronecker's symbol, and the coefficient of the lin- 
ear coupling between adjacent cores is scaled to unity. In 
optics, the discrete equation can be derived by means of 
well-known methods [24j [251 EI] . In the application to 
arrays of plasmonic waveguides, which can be built, for 
example, as a set of metallic nanowires mounted on top of 
a dielectric structure [26] , this equation can be derived in 
the adiabatic approximation, when the exciton field may 



be eliminated in favor of the photonic component (oth- 
erwise, the discrete system will be two-component). It 
is also relevant to mention that the well-known stagger- 
ing transformation [24], u m (t) = (— l) m e~ 2lt where 
the asterisk stands for the complex conjugate, simulta- 
neously reverses the signs of T 2 and B, thus rendering 
the self-focusing and defocusing nonlinearities mutually 
convertible in the discrete system. In particular, the lat- 
ter feature is essential for modeling arrays of plasmonic 
waveguides, where the intrinsic excitonic nonlinearity is 
self-repulsive. In what follows, we fix the signs of Y 2 and 
B by defining T 2 > 0, while B may be positive (self- 
focusing), negative (self-defocusing), or zero. 

As mentioned above, the underlying array can be ac- 
tually manufactured as a uniform one, with all the cores 
doped by an appropriate amplifying material, while the 
HS is singled out by focusing an external pump to a sin- 
gle core. The latter setting is interesting for potential 
applications, as the location of the HS is switchable. 

The model based on Eq. ([2| is the subject of the 
present paper. The Kerr-nonlinearity coefficient, if 
present, may be normalized to B = +1 (self- focusing) 
or B = — 1 (self-defocusing). These two cases are consid- 
ered separately in Sec. \LV\ along with the case of B = 0, 
when the nonlinearity is represented solely by the cubic 
dissipation localized at the HS. 

Dissipative solitons in uniform discrete CGL equations 
were studied by means of numerical methods in Refs. [271 
[28] . We seek analytical solutions for stationary modes 
with real propagation constant k as 



u„ 



ikz 



(3) 



Outside of the HS site, m = 0, the linear lattice gives rise 
to the exact solution with real amplitude A, 



U m = v4exp(-A|ra|), \m\ > 1, 



(4) 



and complex A = Ai +i\ 2 , localized modes corresponding 
to Ai > 0. 

The amplitude at the HS (central) site may be different 
from A, and hence we assume 



= A R + iA T , 



(5) 



for some real constants Ar and Ai. Substituting Eqs. 
([3|, Q and ([5} into the discrete CGL Q yields six non- 
linear algebraic equations for A, Ar, Ai, Ai, A2 and k. 
Straightforward considerations demonstrate that any so- 
lution has Ar = A and Aj = 0, hence the six equations 
reduce to four: 



— 1 + cosh Ai cos A2 = k, 
—7 — sinh Ai sin A2 = 0, 
e" Al sin A 2 - 7 + Ti - EA 2 = 0, 
e" Al cos A 2 - 1 + T 2 + BA 2 = k. 



(6) 



This system of algebraic equations can be solved nu- 
merically by means of the Newton's method for A, Ai, A 2 , 
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perturbed solutions were taken as 




[U m + eV m (z)} e 



ikz 



(7) 
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FIG. 1. Top: The evolution of the pinned-soliton solution ^ 
and Q with A = 0.597, A = 0.626 - 0.846i, and k = -0.202. 
The system's parameters are 7 = 0.5, Ti = 0.9, Y2 = 0.8, 
B — — 1, and E — 0. Bottom: The spectrum of stability 
eigenvalues for small perturbations around this solution. All 
the eigenvalues have non-positive real parts, hence the pinned 
mode is stable. 



and k. For instance, with 7 = 0.5, Ti = 0.9, Y2 = 0.8, 
B = — 1, and E = (these parameters correspond to the 
self-defocusing local nonlinearity and zero cubic loss), a 
physically relevant solution is A = 0.597, Ai = 0.63, 
A2 = —0.85, and k = —0.2. The top panel of Fig. [I] shows 
the stable evolution of the corresponding mode, produced 
by simulating Eq. ([2| with the help of the fourth-order 
Runge-Kutta algorithm, using periodic boundary condi- 
tions. 

Families of pinned modes Q and their stability are 
presented in detail below. Linear gain Ti and cubic 
gain/loss E are used as control parameters in the anal- 
ysis (in the physical system outlined above, their values 
may be adjusted by varying the intensity of the external 
pump). Note that, in the case of B — E = 0, ampli- 
tude A becomes arbitrary and drops out from Eqs. ([6]), 
as Eq. ([2| becomes linear. In this case, T\ may be con- 
sidered as another unknown, determined by the balance 
between the background loss and localized gain in the 
linear system, which implies the structural instability of 
the stationary trapped modes in the linear model. In the 
presence of the nonlinearity, the power balance can be 
adjusted through the value of the amplitude at given Ti, 
therefore solutions, including stable ones, can be found 
in a range of values of the linear gain, Ti. 



III. THE LINEAR-STABILITY ANALYSIS 

The stability of the pinned modes was studied by 
means of the linearization procedure [30 . To this end, 



where V m (z) = X m (z)-\-iY rn (z) is a complex perturbation 
with an infinitesimal amplitude e«l. Substituting this 
into Eq. ([2]), one derives the following linear equations: 



dX n 



1 Y 

2 m 



+£m,o {(T 1 X m - r 2 Y m ) 
-B [2P m Q m X m + (Pi 

-^[(3P^ + g^)x m 4 

dY m 1 „ 



l)Y n 



+ 3Q^) Y m ] 

h IP mQrnYrn\ } 

i)x m 1 



-X 



m+1 



dz 2""^ (/C 
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where P m = Re(C/ m ) and Q m = lm(U m ). An eigenvalue 
problem is obtained by substituting X m = (/> m exp(pz) 
and Ym = ?/; m exp(pz) into Eqs. ([sj) . The pinned mode 
is linearly stable provided that all the eigenvalues have 
Re (p) < 0. An example of the stable numerically cal- 
culated spectrum for the stationary mode considered in 
the previous subsection is shown in the bottom panel of 
Fig. [I} 

IV. NUMERICAL RESULTS 

With the help of the methods outlined above, we con- 
sider three different cases: (i) the self-defocusing non- 
linearity (B = —1), (ii) the self-focusing nonlinearity 
(B — +1), and (hi) zero nonlinearity (B = 0). In each 
case, the cubic gain (E < 0) and loss (E > 0) are inves- 
tigated separately. 



A. The self-defocusing regime (B — —1) 

The top panel in Fig. [2] shows a typical stable soli- 
ton and its eigenvalue spectrum in the self-defocusing 
regime with the linear and cubic gain applied at the 
HS, Ti = 0.9048 and E = -0.16 (we stress that the 
mode is stable in spite of the presence of the unsaturated 
nonlinear gain). With these parameters, Eqs. ([6| yield 
A = 0.6958, Ai = 0.5615, A 2 = 5.2758, and k = -0.3795. 
When the linear gain is increased from F± = 0.9048 to 
Ti = 0.9936, the stable solution attains its largest am- 
plitude, A max = 0.8675, as shown in the middle panel. 
With the further increase of the amplitude, an unstable 
eigenvalue in the spectrum emerges from the origin into 
the right half-plane. The bottom panel depicts such an 
unstable solution with amplitude A = 1.074 > A max , the 
corresponding linear gain being Ti = 0.7731. 

Figure [3] shows amplitude A of the stable (solid) and 
unstable (dashed) pinned modes as a function of linear 
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FIG. 2. (Color online) Examples of stable and unstable 
pinned solitons (left) and the corresponding stability spectra 
(right) in the case of the self-defocusing nonlinearity com- 
bined with the cubic unsaturated gain: B = — 1, E = —0.016, 
7 = 0.5, and Y2 = 0.8. Top: A stable solution found at 
Ti = 0.9048. Middle: The stable solution with the largest 
amplitude, found at Ti = 0.9936. Bottom: An unstable solu- 
tion obtained at I\ = 0.7731 (notice the presence of a positive 
eigenvalue in the spectrum). Here and in other figures, the 
blue solid and red dashed lines represent linearly stable and 
unstable solutions, respectively. 



gain Ti > and cubic gain E < 0. In the absence of 
the cubic dissipation, i.e., at E = 0, there exists a stable 
family of the modes in the region of 0.73 < Ti < 1.11, 
with the amplitude ranging from A = 0.08 to A = 0.89. 
Outside this stability region, any solution governed by 
Eq. ([2| either decays to zero, when the linear gain is too 
small (Ti < 0.73), or blows up to infinity when it is too 
large (T 1 > 1.11). 

Figure [3] shows that, when the cubic gain increases, 
the largest amplitude of the stable soliton and the cor- 
responding linear gain drop, but only by small amounts. 
Naturally, the bifurcation of the zero solution A = 
into the pinned mode takes place at a particular value 
Ti = 0.7286, regardless of the value of E. On the con- 
trary, the unstable branches show a large variation in 
amplitude as E varies. At E = 0, there is a vertical 
asymptote of the unstable branch exactly at Ti = 0.5. 
This implies that, if the local linear gain is too weak, i.e., 
Ti < 0.5, it cannot compensate the background loss with- 
out the contribution from the nonlinear gain. Although 
Fig. [3] shows the unstable branches only in the region of 
Ti > 0, the curves corresponding to E < extend to the 
region of Ti < (linear loss). These unstable solitons 
are supported by the nonlinear gain alone. 

As mentioned in the Introduction, the existence of the 
stability region for the pinned modes in the absence of 
the gain saturation is a remarkable feature of the sys- 
tem. On the other hand, the stability region is (quite 
naturally) much broader in the case of the cubic loss, 



FIG. 3. (Color online) Amplitude A of the pinned soliton as 
a function of linear (Ti) and cubic (E < 0) gain. The other 
parameters are 7 — 0.5, Y2 — 0.8, and B — —1. 

E > 0. Figure [4] shows the respective solution branches 
obtained with the localized self-defocusing nonlinearity. 
When the cubic loss is small, e.g., E = 0.01, there are 
two distinct families of stable modes, representing broad 
small-amplitude (A < 0.89) and narrow large-amplitude 
(A > 2.11) ones. These two stable families are linked by 
an unstable branch with the amplitudes in the interval of 
0.89 < A < 2.11. There is a range of values of the linear 
gain, 0.73 < Ti < 1.13, for which the two stable branches 
coexist, thus making the system bistable. Figure [5] shows 
the coexisting stable solutions in the bistability region. 
In direct simulations, a localized input evolves into either 
of these two stable solutions, depending on the initial am- 
plitude. With the increase of E, the amplitude drops to 
compensate the growing nonlinear loss, stretching the so- 
lution curves in Fig. [4] in the horizontal direction. Simul- 
taneously, the bistability region and the unstable branch 
shrink. Eventually, both of them disappear at E w 0.66. 



B. The self-focusing regime (B — +1) 

For the self-focusing nonlinearity, B = +1, the 
branches of pinned modes ([3| are shown in Fig. [6] as func- 
tions of Ti and E. Under the action of the self-focusing, 
all the pinned modes are unstable without the cubic loss, 
i.e., at E < 0. The left panel of Fig. [7] demonstrates that 
this instability quickly leads to a blowup of the dissipa- 
tive soliton. For the parameters considered here, all the 
unstable solutions originate from point Ti = 0.73. There 
is a vertical asymptote at Ti = 0.5 for the unstable so- 
lutions corresponding to E — 0. These observations are 
identical to those made in the case of the self-defocusing 
nonlinearity (B — —1, see Figs. [3] and [11) 

As shown in the right panel of Fig. [7| the dynamical 
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FIG. 4. (Color online) Solution branches in the case of the 
cubic loss (E > 0). The other parameters are same as those 
in Fig. [3] 
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FIG. 5. The coexistence of stable small- and large- amplitude 
pinned modes (top), and the corresponding evolution of per- 
turbed solutions (bottom) at E = 0.01, in the bistability re- 
gion. Inputs with amplitudes A = 0.3 and A — 2 evolve into 
the small-amplitude and large- amplitude stationary modes, 
respectively. The other parameters are 7 = 0.5, Fi = 1, 
T 2 0.8, and B -1. 



blowup is naturally prevented by the cubic loss {E > 0). 
Figure [6] shows the solution branches for this case too. 
When linear gain Fi falls below a certain threshold, the 
modes do not exist, as the background loss cannot be 
compensated. In this case, all initial condition decay to 
zero. Once Fi exceeds the threshold, the system supports 
the localized modes, which remain stable even at very 
large values of Fi. Figure [8] shows some representative 
examples. For instance, with linear gain Fi = 0.9991, 




FIG. 6. (Color online) Amplitude A of the pinned mode as a 
function of the linear gain (Ti) and cubic loss (E) in the model 
with the self-focusing nonlinearity. The other parameters are 
7 0.5, T 2 0.8, and B +1. 
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FIG. 7. The evolution of perturbed solitons at E = —0.1 
(left) and E = 0.1 (right). The other parameters are 7 = 0.5, 
ri, T 2 = 0.8, and B = 1. 



the system supports a stable pinned mode of amplitude 
A = 2.217 (the top panel). The stable solution with the 
smallest amplitude (A = 0.8521) is found at Fi = 0.6616 
(the middle panel). The solutions with amplitudes A < 
0.8521 are unstable — for instance, the one shown in the 
bottom panel of Fig. [5] 

C. Zero Kerr Nonlinearity (B = 0) 

We have also studied the pinned modes in the case of 
5 = 0, when the nonlinearity at the HS is represented 
solely by the cubic gain or loss, E ^ 0. Figure [9] shows 
the respective solution branches corresponding to differ- 
ent values of E. The linear stability analysis shows that 
all these solutions are unstable in the presence of the cu- 
bic gain [E < 0), while the solutions corresponding to 
the nonlinear loss, E > 0, are always stable. An inter- 
esting feature found here is that the stable and unstable 
branches in the parameter space are symmetric about the 
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FIG. 8. (Color online) Examples of stable and unstable 
pinned solitons (left) and the corresponding stability spectra 
(right) in the case of the self-focusing nonlinearity combined 
with the cubic loss: B = 1, E = 0.1, 7 = 0.5, and T 2 = 0.8. 
Top and middle: Stable solutions found at Ti = 0.9991 and 
Ti = 0.6616, respectively. Bottom: An unstable solution ob- 
tained at Ti = 0.6648. In the latter case, there is a positive 
eigenvalue at p — 0.034. 
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FIG. 9. (Color online) Solution branches as a function of 
the cubic gain or loss E in the case of B = 0. The other 
parameters are the same as those in Fig. [7] 



solution for E = 0. This particular branch has an arbi- 
trary amplitude, as there is no nonlinearity when both B 
and E vanish. All the solutions belonging to this branch 
correspond to linear gain Ti = 0.73, which is again the 
critical value at which the zero background bifurcates 
into the pinned mode (see Figs. [3| [4] and [6|. Several 
stable and unstable modes and their linear spectra are 
shown for the case of B = and different values of E in 

Fig.rroi 



FIG. 10. (Color online) Examples of stable and unstable 
pinned modes (left) and the corresponding stability spectra 
(right) in the absence of the Kerr nonlinearity (B = 0). All 
these solutions have the same amplitudes, A = 3.687. The 
other parameters are 7 = 0.5, and Y2 = 0.8. Top: A stable so- 
lution found at E = 0.02 with Ti = 1.0004. Middle: A stable 
solution found in the linear system, with E — and Fi = 0.73 
(the amplitude is arbitrary in this case). Bottom: An unsta- 
ble solution obtained at E = -0.02 with F ± = 0.4568. In the 
latter case, there is a positive eigenvalue, p = 0.4955. 



D. Stability and instability of the pinned modes 
with respect to finite perturbations 

It has been shown above that, while the pinned modes 
may be stable against infinitesimal perturbations under 
the combined action of the self-defocusing nonlinearity 
(B — —1) and unsaturated nonlinear gain [E < 0), the 
respective stability area is much smaller than that in the 
case of the nonlinear loss (E > 0), see Figs. [3] and 4] 
The apparent fragility of the dissipative solitons in the 
case of E < 0, and their robustness at E > suggest 
an investigation to check the stability of these two types 
of the pinned modes against finite-amplitude perturba- 
tions, which we have carried out by means of system- 
atic simulations of the evolution of the modes perturbed 
by reasonably large initial disturbances where lineariza- 
tion will not be an adequate approximation initial dis- 
turbances. The conclusion is that the "robust" solitons, 
found at E > 0, are completely stable against arbitrary, 
finite amplitude perturbations. On the other hand, one 
can always destroy the "fragile" modes, which are stable 
against infinitesimal perturbations at E < 0, by apply- 
ing perturbations of a sufficiently large amplitude. If, in 
particular, the finite disturbance is applied by suddenly 
making its amplitude larger than it is in the stationary 
solution ("stretching"), the soliton will blow up if the 
stretching factor exceeds a particular critical value. The 
corresponding instability borders for the solitons with 
E = —0.2 and E = are displayed in Fig. [IT 



In addition, the blowup of the linearly stable but "frag- 
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FIG. 11. (Color online) The blue dashed curves show 
the largest initial disturbance, corresponding to the sudden 
"stretch" of the soliton, after the application of which the 
pinned dissipative soliton still relaxes back into the original 
form. A yet stronger stretch quickly initiates a blowup. The 
left and right panels pertain to E — —0.2 and 0, respectively. 
The other parameters being B = — 1, and ri = Y2 = 0.8. 
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FIG. 12. The left panel: The evolution of a strongly per- 
turbed "fragile" pinned mode, at E — —0.2. The right panel: 
the same, but for E — +0.2. Other parameters are as in Fig. 
1 1 1 1 The other parameters are same as those used in Fig. |11| 



ile" pinned mode, caused by arbitrary finite amplitude 
perturbations, and, simultaneously, the full stability of 
the "robust" mode against still stronger perturbations, 
are illustrated by their evolution histories presented in 
Fig. " 
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It may be concluded that the entire space of 
the initial conditions is the attraction basin of the latter 
mode, while for the one supported by the unsaturated 
cubic gain the attraction basin is quite narrow. 



linear site with the linear gain (HS, "hot spot") is em- 
bedded. The system can be readily implemented in the 
form of an array of optical or plasmonic waveguides, ad- 
mitting selective excitation of individual cores, by the 
local application of the pump to the uniformly doped 
cores. Solutions for solitons pinned to the central site 
were found in the implicit analytical form, and their sta- 
bility against infinitesimal and finite perturbations was 
investigated numerically. Stability regions for the soli- 
tons have been identified in the parameter plane of the 
most essential control parameters of the system, viz., the 
linear gain Ti and cubic dissipation E. A nontrivial find- 
ing is a (rather small) stability area for the solitons sup- 
ported by the combination of the local nonlinear unsat- 
urated gain and self-defocusing cubic nonlinearity. On 
the other hand, the combination of the cubic loss and 
self-defocusing nonlinearity gives rise to the bistability 
of the pinned solitons. In the former case, the collapse 
of the linearly stable soliton is caused by finite ampli- 
tude perturbations. These features may be promising 
for potential applications, and call for an experimental 
implementation. 

The work may be naturally extended in different direc- 
tions. In particular, it will be interesting to investigate 
localized modes pinned by pairs of hot spots (cf.the anal- 
ysis of a discrete counterpart of the continuous model of 
Ref. [16 j). A challenging possibility is to develop the 
analysis for modes pinned to the hot spot embedded into 
a two-dimensional linear lossy lattice. In that case, an an- 
alytical solution is not available even for the linear lattice, 
hence the entire analysis should be done in a numerical 
form. 
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